A matched-filter-based reverse-time migration algorithm for ground-penetrating radar data

Ground-penetrating radar (GPR) is a remote sensing technique used to obtain information on subsurface features from data collected over the surface. The process of collecting data may be viewed as mapping from the object space to an image space. Since most GPRs use broad beamwidth antennas, the energy reflected from a buried structure is recorded over a large lateral aperture in the image space. Migration algorithms are used to reconstruct an accurate scattering map by refocusing the recorded scattering events to their true spatial locations through a backpropagation process. The goal of this paper is to present a pair of finite-difference time-domain (FDTD) reverse-time migration algorithms for CPR data processing. Linear inverse scattering theory is used to develop a matched-filter response for the GPR problem. The reverse-time migration algorithms, developed for both bistatic and monostatic antenna configurations, are implemented via FDTD in the object space. Several examples are presented

Radar mapping of Isunnguata Sermia, Greenland

This is the published version. Copyright 2013 International Glaciological SocietyIce thickness estimates using advanced nadir sounding and tomographic radar processing techniques are compared and combined in a study of Isunnguata Sermia glacier, Greenland. Using an ensemble of Operation IceBridge flight lines spaced at 500 m intervals and running approximately along the flow direction, we find there is a statistically excellent comparison between subglacial terrains derived from two-dimensional tomography and gridded nadir sounding. Analysis shows that tomographic data better capture short wavelength (1–2 km) patterns in basal terrain, but interpolated nadir sounding data yield more spatially extensive and continuous coverage across the glacier, especially in deep subglacial troughs. Using derived surface and basal topography maps, we find that driving stress and measured and modeled surface velocity comparisons indicate that basal sliding is an important component of the glacier motion, but is also only weakly coupled to the detailed bed topography save for the deepest troughs. As might be expected for this land-terminating, relatively slow-moving glacier, the subglacial and proglacial topography is similar, suggesting the erosional processes acting on the modern glacier bed once helped sculpt the now exposed land

CReSIS airborne radars and platforms for ice and snow sounding

This work is licensed under a Creative Commons Attribution 4.0 International License.This paper provides an update and overview of the Center for Remote Sensing of Ice Sheets (CReSIS) radars and platforms, including representative results from these systems. CReSIS radar systems operate over a frequency range of 14–38 GHz. Each radar system's specific frequency band is driven by the required depth of signal penetration, measurement resolution, allocated frequency spectra, and antenna operating frequencies (often influenced by aircraft integration). We also highlight recent system advancements and future work, including (1) increasing system bandwidth; (2) miniaturizing radar hardware; and (3) increasing sensitivity. For platform development, we are developing smaller, easier to operate and less expensive unmanned aerial systems. Next-generation platforms will further expand accessibility to scientists with vertical takeoff and landing capabilities

Airborne fine-resolution UHF radar: an approach to the study of englacial reflections, firn compaction and ice attenuation rates

This is the published version. Copyright 2015 International Glaciological SocietyWe have built and operated an ultra-wideband UHF pulsed-chirp radar for measuring firn stratigraphy from airborne platforms over the ice sheets of Greenland and West Antarctica. Our analysis found a wide range of capabilities, including imaging of post firn–ice transition horizons and sounding of shallow glaciers and ice shelves. Imaging of horizons to depths exceeding 600 m was possible in the colder interior regions of the ice sheet, where scattering from the ice surface and inclusions was minimal. The radar's high sensitivity and large dynamic range point to loss tangent variations as the dominant mechanism for these englacial reflective horizons. The radar is capable of mapping interfaces with reflection coefficients as low as –80 dB near the firn–ice transition and as low as –64 dB at depths of 600 m. We found that firn horizon reflectivity strongly mirrored density variance, a result of the near-unity interfacial transmission coefficients. Zones with differing compaction mechanisms were also apparent in the data. We were able to sound many ice shelves and areas of shallow ice. We estimated ice attenuation rates for a few locations, and our attenuation estimates for the Ross Ice Shelf, West Antarctica, appear to agree well with earlier reported results

The trough-system algorithm and its application to spatial modeling of Greenland subglacial topography

This is the published version. Copyright 2014 Herzfeld et al.Dynamic ice-sheet models are used to assess the contribution of mass loss from the Greenland ice sheet to sea-level rise. Mass transfer from ice sheet to ocean is in a large part through outlet glaciers. Bed topography plays an important role in ice dynamics, since the acceleration from the slow-moving inland ice to an ice stream is in many cases caused by the existence of a subglacial trough or trough system. Problems are that most subglacial troughs are features of a scale not resolved in most ice-sheet models and that radar measurements of subglacial topography do not always reach the bottoms of narrow troughs. The trough-system algorithm introduced here employs mathematical morphology and algebraic topology to correctly represent subscale features in a topographic generalization, so the effects of troughs on ice flow are retained in ice-dynamic models. The algorithm is applied to derive a spatial elevation model of Greenland subglacial topography, integrating recently collected radar measurements (CReSIS data) of the Jakobshavn Isbræ, Helheim, Kangerdlussuaq and Petermann glacier regions. The resultant JakHelKanPet digital elevation model has been applied in dynamic ice-sheet modeling and sea-level-rise assessment

Radar sounder evidence of thick, porous sediments in Meridiani Planum and implications for ice‐filled deposits on Mars

Meridiani Planum is one of the most intensely studied regions on Mars, yet little is known about the physical properties of the deposits below those examined by the Opportunity rover. We report the detection of subsurface echoes within the Meridiani Planum deposits from data obtained by the Mars Advanced Radar for Subsurface and Ionospheric Sounding (MARSIS) instrument. The delay time between the surface and subsurface returns is indicative of materials with a real dielectric constant of 3.6 ± 0.6. The real dielectric constant is strongly modulated by bulk density. Newly derived compaction relationships for Mars indicate that the relatively low dielectric constant of the Meridiani Planum deposits is consistent with a thick layer of ice‐free, porous, basaltic sand. The unique physiographic and hydrologic setting of Meridiani Planum may have provided an ideal sediment trap for eolian sands. The relatively low gravity and the cold, dry climate that has dominated Mars for billions of years may have allowed thick eolian sand deposits to remain porous and only weakly indurated. Minimally compacted sedimentary deposits may offer a possible explanation for other nonpolar region units with low apparent bulk dielectric constants

An ultra-wideband, microwave radar for measuring snow thickness on sea ice and mapping near-surface internal layers in polar firn

This is the published version. Copyright 2013 International Glaciological SocietySea ice is generally covered with snow, which can vary in thickness from a few centimeters to >1 m. Snow cover acts as a thermal insulator modulating the heat exchange between the ocean and the atmosphere, and it impacts sea-ice growth rates and overall thickness, a key indicator of climate change in polar regions. Snow depth is required to estimate sea-ice thickness using freeboard measurements made with satellite altimeters. The snow cover also acts as a mechanical load that depresses ice freeboard (snow and ice above sea level). Freeboard depression can result in flooding of the snow/ice interface and the formation of a thick slush layer, particularly in the Antarctic sea-ice cover. The Center for Remote Sensing of Ice Sheets (CReSIS) has developed an ultra-wideband, microwave radar capable of operation on long-endurance aircraft to characterize the thickness of snow over sea ice. The low-power, 100mW signal is swept from 2 to 8GHz allowing the air/snow and snow/ ice interfaces to be mapped with 5 cm range resolution in snow; this is an improvement over the original system that worked from 2 to 6.5 GHz. From 2009 to 2012, CReSIS successfully operated the radar on the NASA P-3B and DC-8 aircraft to collect data on snow-covered sea ice in the Arctic and Antarctic for NASA Operation IceBridge. The radar was found capable of snow depth retrievals ranging from 10cm to >1 m. We also demonstrated that this radar can be used to map near-surface internal layers in polar firn with fine range resolution. Here we describe the instrument design, characteristics and performance of the radar

Bed topography of Jakobshavn Isbrae, Greenland, and Byrd Glacier, Antarctica

This is the published version. Copyright 2015 International Glaciological SocietyThis paper presents the bed topography of Jakobshavn Isbrae, Greenland, and Byrd Glacier, Antarctica, derived from sounding these glaciers with high-sensitivity radars. To understand the processes causing the speed-up and retreat of outlet glaciers, and to enable the development of next-generation ice-sheet models, we need information on bed topography and basal conditions. To this end, we performed measurements with the progressively improved Multichannel Coherent Radar Depth Sounder/Imager (MCoRDS/I). We processed the data from each antenna-array element using synthetic aperture radar algorithms to improve radar sensitivity and reduce along-track surface clutter. We then applied array and image-processing algorithms to extract the weak bed echoes buried in off-vertical scatter (cross-track surface clutter). At Jakobshavn Isbrae, we observed 2.7 km thick ice ∼30 km upstream of the calving front and ∼850 m thick ice at the calving front. We also observed echoes from multiple interfaces near the bed. We applied the MUSIC algorithm to the data to derive the direction of arrival of the signals. This analysis revealed that clutter is dominated by the ice surface at Jakobshavn Isbrae. At Byrd Glacier, we found ∼3.62 km thick ice, as well as a subglacial trench ∼3.05 km below sea level. We used ice thickness information derived from radar data in conjunction with surface elevation data to generate bed maps for these two critical glaciers. The performance of current radars must be improved further by ∼15 dB to fully sound the deepest part of Byrd Glacier. Unmanned aerial systems equipped with radars that can be flown over lines spaced as close as 5 m apart in the cross-track direction to synthesize a two-dimensional aperture would be ideal for collecting fine-resolution data over glaciers like Jakobshavn near their grounding lines

Coherent radar ice thickness measurements over the Greenland ice sheet

This is the published version, also available here: http://dx.doi.org/10.1029/2001JD900183.We developed two 150-MHz coherent radar depth sounders for ice thickness measurements over the Greenland ice sheet. We developed one of these using connectorized components and the other using radio frequency integrated circuits (RFICs). Both systems are designed to use pulse compression techniques and coherent integration to obtain the high sensitivity required to measure the thickness of more than 4 km of cold ice. We used these systems to collect radar data over the interior and margins of the ice sheet and several outlet glaciers. We operated both radar systems on the NASA P-3B aircraft equipped with GPS receivers. Radar data are tagged with GPS-derived location information and are collected in conjunction with laser altimeter measurements. We have reduced all data collected since 1993 and derived ice thickness along all flight lines flown in support of Program for Regional Climate Assessment (PARCA) investigations and the North Greenland Ice Core Project. Radar echograms and derived ice thickness data are placed on a server at the University of Kansas (http://tornado.rsl.ukans.edu/Greenlanddata.htm) for easy access by the scientific community. We obtained good ice thickness information with an accuracy of ±10 m over 90% of the flight lines flown as a part of the PARCA initiative. In this paper we provide a brief description of the system along with samples of data over the interior, along the 2000-m contour line in the south and from a few selected outlet glaciers

Annual Greenland accumulation rates (2009–2012) from airborne snow radar

Contemporary climate warming over the Arctic is accelerating mass loss from the Greenland Ice Sheet through increasing surface melt, emphasizing the need to closely monitor its surface mass balance in order to improve sea-level rise predictions. Snow accumulation is the largest component of the ice sheet's surface mass balance, but in situ observations thereof are inherently sparse and models are difficult to evaluate at large scales. Here, we quantify recent Greenland accumulation rates using ultra-wideband (2–6.5 GHz) airborne snow radar data collected as part of NASA's Operation IceBridge between 2009 and 2012. We use a semiautomated method to trace the observed radiostratigraphy and then derive annual net accumulation rates for 2009–2012. The uncertainty in these radar-derived accumulation rates is on average 14 %. A comparison of the radar-derived accumulation rates and contemporaneous ice cores shows that snow radar captures both the annual and long-term mean accumulation rate accurately. A comparison with outputs from a regional climate model (MAR) shows that this model matches radar-derived accumulation rates in the ice sheet interior but produces higher values over southeastern Greenland. Our results demonstrate that snow radar can efficiently and accurately map patterns of snow accumulation across an ice sheet and that it is valuable for evaluating the accuracy of surface mass balance models